Functional renormalization group for d-wave superconductivity in Hubbard type models 
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The temperature dependence of d-wave superconducting order for two dimensional fermions with d-wave 
attraction is investigated by means of the functional renormalization group with partial bosonization. Below the 
critical temperature T c we find superconductivity, a gap in the electron propagator and a temperature dependent 
anomalous dimension. At T c the renormalized "superfluid density" jumps and the approach to T c from above 
is characterized by essential scaling. These features are characteristic for a phase transition of the Kosterlitz- 
Thouless (KT) type. 



The two-dimensional Hubbard model tlhis often associated 
with high temperature superconductivity [2]. Its ground state 
is assumed to be a d-wave superconductor for an appropriate 
range of doping |3. 4|. Under the hypothesis of d-wave super- 
conductivity we study within a simplified effective model the 
phase transition from the high temperature "symmetric phase" 
to the low temperature superconducting phase by means of 
the functional renormalization group. As characteristic fea- 
tures we find (i) essential scaling as the critical temperature 
T c is approached from above, (ii) a massless Goldstone mode 
for T < T c leading to superfluidity, (iii) a temperature de- 
pendent anomalous dimension r\ for T < T c , (iv) a jump in the 
renormalized superfluid density at T c and (v) a gap in the elec- 
tron propagator for T < T c which vanishes non-analytically as 
T -> T c . 

The features (i)-(iv) are characteristic for a KT phase tran- 
sition [5] which is a natural candidate for a universality class 
with a global U(l) symmetry. Actually, while in the infi- 
nite volume limit the "bare" order parameter vanishes in ac- 
cordance with the Mermin- Wagner theorem, a renormalized 
order parameter [6] spontaneously breaks the U(l) symme- 
try for T < T c , leading to an infinite correlation length for 
the Goldstone mode. Vortices arise as topological defects for 
T < T c . In the Hubbard model the fermionic excitations (sin- 
gle electrons or holes) are an important ingredient for driving 
the transition. Finally, coupling our model to electromagnetic 
fields the photon acquires a mass through the Higgs mech- 
anism while the Goldstone mode disappears from the spec- 
trum, being a gauge degree of freedom. Superfluidity is re- 
placed by superconductivity. Recently, various experimental 
data have been interpreted as signatures of KT phase fluctua- 
tions 171 la, l9l 1 1011 . The direct relation between real high tem- 
perature superconductors and the two-dimensional Hubbard 
model is difficult, however. 

The problem of d-wave superconductivity in the Hubbard 
model can be separated into two qualitative steps in the renor- 
malization flow. In the first step the fluctuations with momenta 
IP 2 ~~ I > A 2 (with pf on the Fermi-surface) have to gener- 
ate a strong effective interaction in the d-wave channel. This 
phenomenon has indeed been found in renormalization group 
studies (IlIIlEHISEISBIIl, the d-wave cou- 
pling being triggered by spin wave fluctuations. It also has 
been found in various other approaches, see e.g. I20I, I2TI l22tl . 
The second step, involving momenta |p 2 — p^.| < A 2 for the 
fermions and p 2 < A 2 for bosonic bound states, has to deal 
with a situation where the d-wave channel coupling domi- 



nates. We only address here the second step. We find that 
the universal behavior in the vicinity of the critical temper- 
ature is dominated by effective bosonic fluctuations and be- 
comes independent of the microscopic details. This justifies 
our approach with a single coupling Aj at the scale A in the 
d-wave channel, even though this remains, of course, only an 
approximation. The prize to pay is that A4 is only an effec- 
tive coupling - its relation to the parameters t and U of the 
Hubbard model is expected to depend on T and the chemical 
potential fi. 

We emphasize that we obtain a KT-like phase transition 
within a purely fermionic model. The effective bosonic 
degrees of freedom arise as composite fields - we do not 
start from a bosonic effective theory. The transition from a 
fermionic to a bosonic description is a result of the renor- 
malization flow. The complete decoupling of the fermionic 
degrees of freedom for the description of the critical behav- 
ior is a non-trivial result which only holds for a sufficiently 
large critical temperature T c . (For T c — > interesting modi- 
fications reflect the different universality class of a quantum 
phase transition.) Several new features are directly related 
to the dominance of fermionic fluctuations sufficiently away 
from the critical region: We find that the critical temperature 
T c is below a "pseudo-critical temperature" T pc where the ef- 
fective interaction strength for the electrons diverges. On the 
other hand, the temperature range T c < T < T pc is strongly af- 
fected by fluctuations of composite bosons. For the transition 
to d-wave superconductivity we find that T c is actually rather 
near to T pc . This contrasts to the transition to antiferromag- 
netism foru = where T c has been found substantially below 



T pc |23|24H. Other fermionic aspects concern the relation be- 
tween the gap in the electron propagator and the superfluid 
order. We also find a quantitatively important finite size effect 
which modifies the essential scaling for T — > T c . Again it is 
dominated by the fermionic fluctuations. 

Our starting point is an action with nearest and next-nearest 
neighbor hopping on a two dimensional square lattice with an 
attractive interaction in the d-wave channel 

T*W = ^HQXiaiF + ^-fxMQ) (1) 

Q 

- A d 2_j 8 Ki+K2 , K3+Ki d((\ii-\i2)/2) 

K\,..,,K4 

We employ the short-hand notation Q = (a>f, q), top = (2n + 



2 



T E« j" n -firtp- The Fermi surface is spec- 



if, n e Z, £ e - . ^j.^ 
ified by e q = -2f(cosgi + cos ^2) - 4?' cosgi cos ^2 = and 
of(q) = cos q\ - cos q2 is a d-wave form factor. All lengths are 
measured in units of the lattice distance which we set to unity. 
We use partial bosonization with a bosonic field (p standing 
for a fermion bilinear ~ ipi/t in the d-wave channel. For the in- 
tegration of the fluctuations we employ the concept of the ef- 
fective average action H25I1 . which amounts to the quantum 
effective action (generating lPI-vertices) in presence of an in- 
frared cutoff k. It interpolates between the microscopic action 
r^co ~ S and the full quantum effective action = T. 

The flow of the effective av erag e action obeys an exact 
renormalization group equation 112611 . 
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0*1*0, ijt] = ^STr {d k In [Tf\</>, ift] + R k ]} , (2) 

where Rk denotes an infrared cutoff and dk = 2,-(3*RiL) «tr (the 

index i refers to cutoffs for fermions and bosons). The "super- 
trace" STr runs over field type, momentum and internal in- 
dices, and has an additional minus sign for fermionic entries. 
The functional differential equation (f2]i involves the logarithm 
of the full inverse propagator FJ. [<f>, ift] in presence of back- 
ground fields (second functional derivative of T^), regulated 
by Rk- Approximations to the solution of eq. (O proceed by a 
truncation of Ta on the r.h.s. with a suitable ansatz. 
Our truncation for the effective average action is 

Ft = J]iA + (0(/w F + e q -^(0 (3) 

Q 

+ Yj f(QK iZc °B + A [q] 2f2 ¥(G) + X v ffl 

Q x 

K,Q,Q' 

x(<f WtCGWGO - W^(0^(6')) . 

r 1 

with X - (t, x), YjX - J T dr and for bosons Q - {u>b, q), 
o>b = 2mnT, m e Z. Further [q] 2 is defined as [q] 2 - q\ + q\ 
for qi e [-Ji, n] and continued periodically otherwise. We 
approximate the effective potential by a quartic polynomial 



U((f>) 



m 2 5+^6 2 (SYM) 
~%(6-6 Q ) 2 (SSB) 



(4) 



with 6 = <p*<f>. Here (SYM) is used as long as the minimum 
of U is located at <f> = 0, while the regime with spontaneous 
symmetry breaking (SSB) is characterized by non vanishing 
60 = 4>Q<f>o- The order parameter <f>o determines the size of 
the gap in the fermion propagator A(q) = \d(q)<f>o\. The gen- 
eralized couplings Z, A, ill 2 , 60 and 1^ depend on k. At the 
scale A our truncation is equivalent to the fermionic action ([TJi 
provided 



(m/?) 2 = -^-, 1^ = 0, Z = A=0, 
as can be seen by solving for <f> as a functional of ift. 



(5) 



We will follow the flow until a physical scale k p h where 
k~l corresponds to the macroscopic size of the experimental 

probe. (In practice, our choice k p k/t = 10~ 9 = e~ 20 - 723 corre- 
sponds to a probe size of roughly 1 cm.) For T < T c we find 
a non vanishing superfluid density 6o(k p i,) and therefore a non 
vanishing gap A, while So(k P h — > 0) = 0. 

In addition to the truncation we have to specify the regulator 
functions Rk- For the fermion fluctuations we use a tempera- 
ture like regulator, 



<(G) = *W-(^-l), with T 4 k =T 4 +k 4 



(6) 



For k > T the cutoff R^ in the inverse fermion propagator 
suppresses the contribution of all fluctuations with momenta 
|p - p/H 2 < (nk) 2 , even for T = 0. It becomes ineffective for 
k <K T where no cutoff is needed anymore. The bosonic cutoff 
regularizes the fluctuations of long range bosonic modes. We 
use a "linear cutoff" for the space like momenta j27il28ll 



R B k (Q) =A-(k 2 - [q]¥) ©(k 2 - [q]¥) . 



(7) 



The flow equations for the couplings are obtained by pro- 
jection from ©. In (SYM) the effective potential flows as 



d k Uk(6) 



2T(d k T k ) 



r — 



-y a tanhy t5 



(8) 



s 



d k Rt(Q)[Pf{-Q)+m 2 



Q L * 
where 



[Pf (0 + in 2 + 2^<5][Pf(-0 + in 2 + 2A^8] - A 2 6 2 ' 



= ^L^-^+A^q), 



Pf(Q) = iZaj B +A[q] 2 t 2 +R B k (Q). 



7S 



(9) 



For the mass parameter and the bosonic quartic coupling this 
yields 



- 2 " 
do 



k*4 



°1_ 
d 2 5 



d k Uk(6)\ s= o. (10) 



The fermion fluctuations generate kinetic and gradient 
terms for the bound state bosons (SYM regime) 



d k Z = -T(d k T k ) 



r d 2 q , d d 
J_„ (2ny dT k wo) 

tanh^ +tanh(^ + §) 



IT' 



and 



1 d 

d k A = -27WjfcW 1 



a 



t 2 dT k d(P) J_ 



f— 



^{q - l -^) 2 



X- 



tanh(^) + tanh(^-^) 



Tk(e q + e/ ei -q - 2fi) 



j=0 



(12) 



There are no bosonic contributions to the flow of Z and A to 
this level of approximation. The flow starts in the symmetric 
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regime © with initial values ((5]) and we choose A/f = 0.5. 
For temperatures below a pseudo-critical temperature T pc the 
scale dependent mass term m 2 vanishes at a certain critical 
scale kssB- This signals local ordering and indeed corresponds 
to a divergent effective four fermion coupling Ad(k). Below 
kssB we use the truncation (SSB) in ®. 

We do not consider here the region of very small T where 
interesting quantum critical phenomena may occur [ 29] . Then 
classical statistics dominates the flow for k <K T. A typi- 
cal flow in the SSB regime is depicted in fig. [TJ where we 
show k = t 2 A6o/T for various values of T and t' jt - -0.1, 
Ad/t = 5/8, fi/t = -0.6. The fermion fluctuations (first term 
in eq. l[8]l) first drive m 2 (k) to zero and induce non vanishing 
K(k), starting from K(k SSB ) = 0. For k <K kssB the effect of the 
fermion fluctuations is strongly reduced due to the suppres- 
sion with powers k/T from d k T k and the presence of a non 
vanishing gap A. Now the bosonic fluctuations take over and 
tend to reduce the value of K{k) as k is lowered further. For 
the region around and on the right of the maximum of K(k) 
in fig. [TJwe have already k <K T. Then the Matsubara fre- 
quency m — dominates the bosonic contributions and the 
system gets reduced to classical statistics for the (9(2) linear 
cr-model. One can see this effective reduction from 2 + 1 to 2 
dimensions explicitly by evaluating the bosonic contribution 
in eq. ([8]) (the second term, denoted by d k Uf). Keeping only 
the m = contribution the flow equation simplifies consider- 
ably 



„ TAk r" d z q , q z t z , , 

8 k U B k {5)=— J _^(2-^)(l-i_))0(£ 2 -q 2 f 2 ) 

(13) 



— + 



Ak 2 + (S- S )Xf, Ak 2 + (36 - 6 )A 4 , 



Here rj(k) = — defines the anomalous dimension. 

For a better understanding of the scaling behavior we intro- 
duce rescaled and renormalized quantities 



Tk 2 



U . 



t 2 A 



t 2 A 



(14) 



This yields for the flow of the rescaled effective potential, now 
at fixed 5 instead of fixed 6, 

~ du k t 2 

kd k u k = -2u k + rj(k)6—=- + —Tkd k U k \ s , (15) 
86 i k 

implying for the renormalized field expectation value k and 
bosonic quartic coupling 

1 d d 2 
dkK = - — -~d k u k \ , dkAj, = —rd k u k \ . (16) 
A$ do o z o 

The bosonic contribution to the anomalous dimension in SSB 
is the same as for the linear (9(2) model in two dimensions 



1 AiK 
r,(k) = ~- 



7T (1 +2,V) 2 



(17) 



while the fermionic contribution can be neglected. 

The /?-functions can be expanded in powers of (IA^kY 1 . In 
particular, to leading order one finds 

^2 

kd k K =p K = 0(k- x ) , kd k A = -lA^ + -± + 0(k- x ), (18) 

resulting in an infrared attractive fixed point for A$. The run- 
ning of k can be mapped to the running of the cou pling i n 
the non-linear cr-model for arbitrary (9(A0-models [6, 30, 31]. 
For A = 2 also the term ~ k~ 1 in /3 K has to vanish, even though 
this is not seen in the truncation (O. The flow of the classical 
linear (9(2)-model has already been studied with much more 
extended truncations [31], taking into account an arbitrary in- 
dependence of u(5) and arbitrary 5-dependent wave function 
renormalizations, e.g. A(6). All these functions flow quickly 
to a scaling form which depends only on one parameter k. 
Evaluated for the scaling solution one finds to a very good 
approximation for k > kkt 



kd k K — f3 K — a\K» - k\ 2 &{k, - k) 



(19) 



with k„ — 0.248 and a — 2.54 m. 

We will use this improved 
truncation for f3 K for k < 77 12 for the region k > kkt =0.1. 
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FIG. 3: Temperature dependence of d-wave superconducting order. 
We plot the renormalized superfluid density k and the fermionic gap 
A(0, |)/f. Solid/dashed lines are based on eq. lfT9ll/ll3l. 



This results in the solid lines in fig. [T| whereas the truncation 
© gives the dashed lines. 

We define the critical temperature T c as the temperature for 
which the field expectation value a- just vanishes at the physi- 
cal scale kph/t = 1CT 9 , see fig. [1] For T < T c one finds a su- 
perfluid condensate in a probe of macroscopic size of roughly 
1 cm. The dependence of T c on the specific choice of k p % is 
extremely weak. For our parameters t = 0.3 eV corresponds 
to T c = 128 K while for alternative Aj/t = 1/2, 1 we find 
T c = 58.3, 366 K, respectively. We also note that T c is only 
mildly affected by the improvement of the truncation ( fT9b . 

A specific feature of the KT phase transition is essen- 
tial scaling for the temperature dependence of the correlation 
length £ = m R l = (m z /A)~2 just above T c , 



m R = ce VW-->-< . (20) 

The finite size correction C(k p h) ~ (1 + §k« 2 In j— )~ 2 vanishes 
logarithmically for k p f, — > 0. This follows from eq. ( fT9l ). 
Close to the critical temperature the omission of the finite size 
correction (i.e. C(k p h = 0) yields to a substantial derivation 
from the straight line in fig. 2. 

For T > T c one finds that K(k) reaches zero at a scale ksR > 
k p h. Continuing the flow in the symmetric regime for ksR > 
k > kph yields m« = 0.50 ksn- We show iriR(T) in fig. [2] 
(same f', Ad, [i in all figs.) and find that essential scaling is 
smoothened due to k p i, > 0. 

From fig. 1 we can also extrapolate the pseudo-critical tem- 
perature T pc . This corresponds to the temperature where k{K) 
bends over immediately after becoming non-zero at k = kssB 
and is driven to zero again more or less at the same scale 
ksR ~ kssB- In other words, T pc corresponds to the high- 
est temperature for which the initial flow in the symmetric 
regime hits a vanishing mass term m 2 and therefore a diver- 
gent effective four fermion vertex ~ ;?T 2 . Often this tem- 
perature is associated with the critical temperature. We find 
Tpc/t - 0.0372, only slightly above the value of the critical 
temperature T c /t = 0.0368. 




0.7 0.75 0.8 0.85 0.9 0.95 1 

T/T c 

FIG. 4: Temperature dependence of the critical exponent rj. The solid 
line is based on eq. l |19t , the dashed line on eq. J3J. 

It is characteristic for two dimensions that the renormal- 
ized order parameter k jumps at T c even though the transition 
shows scaling behavior. (Strictly speaking, this discontinuity 
only occurs for k p i, — > whereas it is smoothened for nonzero 
k p h.) The scaling solution at T c corresponds to a fixed point k, 
where /3 k (k,) = with /3 k (k < /c«) > 0. In the region K(k) < k* 
the flow drives k to zero and the system is in the symmetric 
phase. In contrast, for k > k* the flow can never cross the 
fixed point. Therefor K(k p h) remains larger than k*, resulting 
in a minimal value for k in the ordered phase. For k p h — * this 
results in a jump Ak = /e* between the two phases. In super- 
fluids K(kph) determines the renormalized superfluid density. 
The value resulting in the vortex picture [5] is k* = -. Our 
truncation ( fl9l ) is already parameterized in terms of k„. In fig. 
[3] we plot the value of the superfluid density K(k p h) as well as 
the gap in the electron propagator A(0, |) in a fixed momen- 
tum direction. In the present truncation A would vanish for 
kph — * 0, whereas k would remain nonzero. In order to settle 
this issue more precisely one would have to follow the flow of 
the Yukawa coupling that we neglected here. For t = 0.3 eV 
and T = T c /2 we obtain A(0, §) = 7.7 meV. 

We finally turn to the anomalous dimension for T < T c . For 
weakly inhomogeneous fields (low q 2 ) the gradient terms in 
the effective theory (Landau theory) for the composite bosons 
(f> take the form 

r = j ^f(q)A(q 2 , kph) q 2 c/,(q) . (21) 

Our truncation ^ has not yet taken the q 2 -dependence of A 
into account: the flow equations apply for A(k) = A(q 2 = 
0, k). Nevertheless, the momentum dependence of A can be 
easily inferred by noting that q 2 also acts as an infrared cutoff 
such that Rit becomes ineffective for k 2 < q 2 . A reasonable 
approximation uses SjA(q 2 , k) = (5jA(q 2 = 0, k)&(k 2 - q 2 ). 
The solution at k p h reads for q 2 > k 2 ph , q = ^[^, 

plnq 

A(q 2 ) =A(q/q)-'> , r] = d In k i](k) / ln(q / q) . (22) 

J\n q 

In fig. [4|we show 77 for qlq-e 5 and observe the characteris- 
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tic dependence on T. As a result the static correlation function 
(<p{r)(p*(0)) c decreases ~ f~ n for large r. For a superconduc- 
tor one replaces in eq. d2TT i q 2 by -(V - 2z'|e|A) 2 with A the 
space components of the electromagnetic potential. Inserting 
4> = <po this results in an anomalous Landau theory for the 
magnetic field B = d\A2 — diAi. Instead of the photon mass 



term ~ A 2 responsible for ordinary superconductivity in three 
dimensions we now find a term |A| 2 ~ n . This may open a win- 
dow for a measurement of the anomalous dimension. 
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